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Abstract 

The appearance of sharp vorticity gradients in two-dimensional hydrodynamic turbulence and 
their influence on the turbulent spectra is considered. We have developed the analog of the vortex 
line representation as a transformation to the curvilinear system of coordinates moving together 
with the di-vorticity lines. Compressibility of this mapping can be considered as the main reason 
for the formation of the vorticity discontinuities at high Reynolds numbers. For two-dimensional 
turbulence in the case of strong anisotropy the vorticity discontinuities can generate spectra with 
the fall-off at large k proportional to A;^'^ resembling the Kraichnan spectrum for the enstrophy cas- 
cade. For turbulence with weak anisotropy the k dependence of the spectrum due to discontinuities 
coincides with that of the Saffman spectrum: k~'^. We have compared the analytical predictions 
with direct numerical solutions of the two-dimensional Euler equation for decaying turbulence. We 
observe that the di-vorticity is reaching very high values and is distributed locally in space along 
piecewise straight lines. Thus, indicating strong anisotropy and accordingly we found a spectrum 
close to the /c~^-spectrum. 
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I. INTRODUCTION 



This paper is concerned with investigations of two-dimensional (2D) hydrodynamical tur- 
bulent flows. In particular, we study the formation and dynamics of very sharp vorticity 
gradients and their influence on the energy spectrum in the enstrophy cascade regime. We 
may consider two kind of turbulent spectra. The first one was suggested by Kraichnan in 
1967 1], it corresponds to the enstrophy cascade directed to the small-scale region where 
viscous dissipation becomes essential. The Kraichnan spectrum follows, up to the logarith- 
mic factor (see P]), a power law for the scales intermediate between source and sink (the 
inertial interval): E(k) ~ rf/'^k^^ where rj is the enstrophy dissipation rate. Recall that 2D 
turbulence additionally is characterized by an inverse energy cascade for large-scales leading 
to E{k) ~ k~^f'^ (see e.g. [1]). However, in the present paper we will only be concerned 
with the small-scale region of the spectrum. The second spectrum suggested by Saffman in 
1971 PI yields another power dependence: E{k) ~ k~^. According to Saffman, in decay- 
ing 2D turbulence vorticity discontinuities (in absence of viscosity) will form because fluid 
elements with different values of vorticity will be driven close together by the flow. Due 
to vorticity conservation the appearance of discontinuities will provide the conservation of 
all other invariants involving vorticity, / u'^dS, n = 3,4, .... Accounting for a finite viscos- 
ity Saffman considers the "discontinuities" to have a small width 6, which results from the 
balance between inertial and viscous forces. At high-Reynolds number this size is assumed 
to be very small in comparison with the length along the discontinuities, L, which may be 
assumed to be of the same order as the characteristic energy-containing length-scale. Under 
the assumption of isotropy and a dilute distribution of discontinuities Saffman suggested 
that the energy spectrum at large k could be constructed as a superposition of the spectra 
from the individual discontinuities resulting in: E{k) ~ k~^ . 

From the first sight, the spectra obtained by Kraichnan and Saffman look like two different 
answers, but indeed, as we show in this paper, it is possible to establish some connection 
between them. This may be seen from the Fourier transform of a step function. Let us 
assume that the vorticity has a jump F = F(x) along the line y = and at first neglect 
effects connected with bending of the line. Then we can write: 

Hence it is immediately seen that the Fourier transform will have a power-law fall-off at large 



k, i.e., inversely proportional to ky multiplied by some function of k^ due to dependence 
of r on X. If we neglect the dependence on kx replacing it by some constant, then we 
immediately obtain an energy spectrum with a power dependence similar to the Kraichnan 
spectrum: E{k) ~ k~^. This is an important conjecture demonstrating that a spectrum 
similar to the Kraichnan spectrum, which is often observed in high resolution numerical 
simulations, may be related to discontinuities of vorticity which can be considered as possible 
candidates for singularities in ideal fluids in two dimensions. However, for viscous fluids, 
i.e., in the framework of the Navier-Stokes equation, such singularities are impossible: the 
initial smoothness of the solution will remain as proven first by Olga Ladyzhenskaya many 
years ago 0]. Within the 2D Euler equations for incompressible fluids the vorticity is a 
Lagrangian invariant and can never be singular, but its gradient might, in principle, become 
infinite in a finite time. Up to now this is an open question. 

It is necessary to mention some examples of 2D flows presented by Yudovich j^, where the 
appearance of weaker singularities (vorticity is allowed to be discontinuous but bounded) are 
possible, however, they are formed in infinite time. Another approach based on the numerical 
analysis of the complex singularities for the inviscid flow with two-mode initial conditions 



showed that the width of its analyticity strip follows a ln(l/t) law a t s 
Additionally, many numerical experiments for 2D turbulence (see. 



lort times 
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) show that with 



a good accuracy the Saffman spectrum is formed at the initial stage, before the excitation 
of the long-scale coherent vortices. The high-resolution numerical simulation performed 
by Legras et al. [12 demonstrated the power dependence k~^'^. Analytical calculations 
presented by Gilbert 1^ using arguments based on the existence of spiral structures give a 
power dependence with exponent between —3 and —4 (see also ll^ and In particular, 

we would like to point to the very interesting paper by Ohkitani (l6|, where by means of the 
Weiss decomposition 3| it was shown that the so-called /i-regions {h- hyperbolic in the sense 
of Ref. i-e., regions where straining is dominating over vorticity) give the spectrum k~^, 
i.e., coinciding up to a logarithmic factor with the Kraichnan spectrum, the contribution from 
the e-regions (e- elliptic, i.e., vorticity dominated regions) yield the Saffman s pec trum ~ 
Note, that similar ideas based on the wavelet analysis were developed in ^ to separate 
2D turbulent flow into regions having different dynamical behaviors. The appearance of a 
power type spectrum in the short-wave region has been connected with different physical 
mechanisms like vortex merging j2^, [2^ and vortex stripping [2^, [2^ which give a certain 
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confirmation of tlie original idea of Saffman P]. 

In this paper we present some qualitative physical arguments in favor of the formation 
of vorticity discontinuities in the 2D Euler equations for smooth initial conditions. The 
main idea in the description of 2D flows is connected with using the vorticity as Lagrangian 
invariants. Kuznetsov and Ruban |24| (see also |25|]) developed a new kind of description 
for three-dimensional vortical flows - the so-called vortex line representation (VLR). This 
representation is based on the mixed Lagrangian-Eulerian description and connected with 
movable vortex lines. The VLR, which is a mapping to a curvilinear system of coordinates, 
turns out to be compressible, this is considered to be the main reason for breaking in 
hydrodynamics. Here we demonstrate how this approach can be modified for 2D Euler 
hydrodynamics. The main obs_ervation is that for 2D flows the curl of the vorticity, sometimes 

t satisfies the same 



referred to as the di- vorticity |9 



equation as, e.g., the equation for the magnetic field in ideal 2D magneto- hydrodynamics 
(MHD). Therefore the generalization to the 2D Euler equations becomes straightforward. In 
the local case, as it was demonstrated for 2D MHD in Ref. the vorticity plays the role 
of a Lagrangian coordinate and the other variable coincides with the Cartesian coordinate, 
say, X. In terms of these variables 



j the 2D Euler equations transform into equations 
of motion for a layered fluid, similar to stratified fluid, where each layer is labeled by its 
vorticity u. In terms of the new variables the "new" hydrodynamics becomes compressible. 
The derivative y^j plays the role of density of each layer, as a function of time and coordinate 
X. This characteristics is proportional to the width between two neighboring layers with 
closed vorticity contours. 

Another aim of this paper is to revisit the energy spectra for 2D turbulence with empha- 
sis on the angle distribution, following the arguments of Saffman connected with vorticity 
discontinuities. Using the stationary phase method we demonstrate that the contribution 
from one discontinuity is very anisotropic: it has a sharp angular peak along the direction 
perpendicular to the discontinuity. In the peak the energy spectrum falls-off like at large 
k. After average over angles in the case of isotropic turbulence the spectrum coincides with 
the Saffman spectrum 

In order to check whether the spectrum at large k is defined by vorticity discontinuities we 
have performed numerical experiments on decaying turbulence based on a direct numerical 
solution of the 2D Euler equations. In the turbulent state when the formation of power tails 
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is observed we examine the structure of the di-vorticity. We found that the di-vorticity is 
distributed very sharply in space concentrated on a random net of hues. In our opinion, 
these results can be interpreted in favor of the Saffman mechanism for the formation of 2D 
turbulent spectra due to discontinuities. 

II. TWO-DIMENSIONAL ANALOG OF THE VLR 

Consider a 2D ideal fluid, described by the Euler equation for the vorticity Lj{x,y,t), 

|^ + (v-VV = 0, divv = 0, (1) 
where the velocity field v defines the vorticity: 

a; = V X V 



dx dy 

Equation (^Q) shows that the vorticity is a Lagrangian invariant advected by the fluid, i.e., 

u!{x, y, t) = inv 

along a fluid particle trajectory defined as solution of the system of ordinary differential 
equations (ODE's), 

— = v(r,t), r\t=o = a. (2) 
Let us introduce the divergence-free vector field B with the components 

dy' ^ (9x ' 

i.e., B = curlcui. It is easily to see that this vector is tangent to the line uj{x,y) = const 
because the vorticity gradient Vcu = {d^uj, dyu) is normal to this line. The equation of 
motion for the vector B can be obtained from Eq. after differentiating with respect to 
coordinates: 

9B 

— = curl[vxB]. (3) 

Thus, the vector B constitutes a frozen-in quantity. Sometimes, it is called as the di-vorticity 
(see jsl). By introducing the substantial (material) derivative, d/dt = d/dt + (v ■ V), Eq. 
(13)) can be rewritten as 

^^(B.V)v. (4) 



Hence, we observe that |B| will locally increase due to stretching of the di-vorticity lines, 
i.e., when 

-— = (B.SB)>0. (6) 

where 

^ 1 / dvk dvi \ 
2 \dxi dxk J 

is the stress tensor. Increasing (or decreasing) the di-vorticity field, based on the equation 
dni), is not sufficient to clarify the physical mechanism for its growth. As is seen from 
Eq. only one velocity component, v„, normal to the vector B changes the field B. In 
this case the tangential component (parallel to B) plays a passive role providing the 
incompressibility condition: divv,- + divv„ = 0. This observation is the key point for 
introducing the vortex line representation (VLR) for the three-dimensional Euler equations 
(see, e.g. |2^). To construct the analog of VLR for the 2D Euler equations we consider new 
Lagrangian trajectories, given by the v^, 

(It 

— = Vn{r,t), r\t=o = a. (6) 

The solution of these ODE's defines a new mapping 

r = r(a,t) (7) 

which is different from that given by Eq. In terms of this mapping the di-vorticity 

equation ^ can be integrated (for details see, e.g. 

B(M) = (8) 

where Bo(a) is the initial di-vorticity, J is the Jacobian of the mapping (|7j): 

J ^ d{x,y) 

According to the definition of this mapping its Jacobian is not fixed, it may change in time 
and space. In other words, the mapping r = r(a, t), as a change of variables, represents a 
compressible mapping. This means that the di-vorticity lines can be compressed. In this 
approach the velocity of motion of di-vorticity lines is nothing else than the normal velocity 
v„. 

It is interesting to note that this approach in slightly different form what was suggested 
in j3| • In this paper the basis of the approach is the mixed Lagrangian- Eulerian description 



when all desired quantities are considered as functions of vorticity uo (or any other Lagrangian 
invariant) and a Cartesian coordinate x. 

The VLR given by (jHI), (0), (jSj) with the local change of variables r = r(a, t), does not 
work at singular points where the B-field vanishes, 

B(r(t),t) = 0. (9) 

and where, respectively, the normal velocity is not defined. For vorticity uo these points are 
nothing more than maximal, minimal or saddle points. It is easy to see that the null points 
are advected by the fluid, but the velocity v at these points is defined through the B-field by 
inverting the Laplacian operator: v = — A^^B. The null-points for the normal vector field 
n(r) represent topological singularities. Topological constraints as additional conditions to 
the system (0), (jH}, are written as integrals of the vector field n(r) along a loop enclosing 
the null-points: 

j{V^ ■ dr) = 27rm, (10) 

where cp is the angle between the vector n and the x-axis and m, being a topological charge, 
is an integer equal to the total number of turns of the vector n while passing around the 



closed contour with the null-point inside it (see also 28|)- For instance, for X-points or 
0-points, m = ±1 . 

As well known from our knowledge in gas-dynamics compressibility of the mapping is 
a main cause for steepening and ultimately breaking, resulting in the formation of sharp 
gradients for the velocity and density of the gases. This happens in finite time and in the 
general situation the singularity first appears in one separate point, i.e., it may be related 
to collapse. In gas-dynamics this process is completely characterized by the mapping deter- 
mined by the transition from the Eulerian to the Lagrangian description. Vanishing of the 
Jacobian corresponds to the emergence of a singularity. For three-dimensional incompress- 
ible Euler equations compressibility of the VLR is a possible reason for appearance of infinite 
vorticity in one separate point that results in breaking of vortex lines. The first study of 
vortex-line breaking for three-dimensional integrable hydrodynamics with the Hamiltonian 
Lj\dr was performed by Kuznetsov and Ruban |23]. Recent numerical experiments 



have confirmed the possibility of this type of scenario. 
The Jacobian in dominator of the expressions (jH)) can become zero, which will result in 
infinite value of the di-vorticity. We do not see any restrictions by which this process can 



be forbidden. In 2D hydrodynamics, however, compressibihty of the mapping guarantees 
only compression of di-vorticity hnes corresponding to the formation of sharp gradients for 
vorticity. Probably, the breaking process in 2D happens in infinite time (see, for instance, 
P]). The most important point for us is that the tendency indeed does exist and it is possible 
to imagine that this process may be accelerated in the presence of external forces driving 
the turbulence. 



III. 2D SPECTRUM 



In the previous section we gave some arguments in favor of formation of sharp gradients 
of the vorticity in 2D Euler flows. Everywhere below we will suppose that this process is 
possible and consider how it can effect the form of turbulent spectrum. For 2D turbulence, 
in the presence of finite viscosity and external forces, we will assume that the sharp vorticity 
gradients have a finite value inversely proportional to the characteristic width of discontinuity 
6, which is defined from the balance between inertial and viscous terms. At high Reynolds 
number the width 6 will be much less than the characteristic (energy-containing) scale L. 
In the turbulent state such discontinuities are naturally assumed to form a set of vorticity 
gaps with random positions of their centers, random forms and random distributions over 
angles. Our aim is to calculate the contribution to the spectrum from such discontinuities. 
We will be interested in the region of wave numbers k lying between and the inverse 
width 6'^: 

To simplify the problem all gaps are supposed to be concentrated on pieces of straight lines 
(finite intervals) with vorticity gaps vanishing at the endpoints of the intervals. As we will 
see later this simplification is not so essential. The answer, which we will get, will also 
account for bends of the discontinuity lines. 

To find spectrum we need to calculate the Fourier transform from of pair correlation 
function: 

F(r) = (cu(x)cu(x + r)), 
where angle brackets means average over the ensemble of discontinuities. Hence the energy 
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density spectrum e(k) is given by the standard formula: 




where uik is the Fourier transform of the vorticity <^(r), 



the over-bar denotes average with respect to random variables, and S is the average area, 
which is assumed to be sufficiently large. 

Consider first one discontinuity with the center at Yq — (xq, yo) oriented along the y-axis. 
Then for the y-derivative of uj we have. 



Here r(,T) is a continuous function of x inside the interval [xi, X2\ vanishing at the endpoints 
X = Xi^2 and equal zero outside the interval. 

Hence, the Fourier transform from the singular part of ou is given by the integral: 



where k = {k^, ky). This is the contribution from one discontinuity. If we assume that the 
discontinuities are not very densely distributed, they may be considered "independent" and 
the spectrum for the whole ensemble of discontinuities may be obtained by a superposition 
of the spectra from the individual discontinuities, i.e., from the summation with respect to 
all discontinuities which results in 



Here we have introduced two unit vectors: normal Hq, and tangent Tq {na^ — — 
1, {jO-oJa) = 0) characterizing the orientation of the a-th discontinuity. The coordinates 
a^ia, a;2a, I/oq together with the two unit vectors define completely the positions of the dis- 
continuities. 

To find an enstrophy spectrum one needs to perform average of \lo\}\' with respect to 
all random variables. Assuming the coordinates of the discontinuities to be randomly dis- 
tributed uniformly in space, the first average gives: 



dy 



T{x)8{y — yo) + regular terms. 



(11) 





(12) 




(13) 
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Here is the number of discontinuities in area S, angle brackets means the average with 
respect to both xi,X2 and angle distribution. 

Since we are interested in short-wave asymptotics of the spectrum, kL ^ 1, the integrand 
in (jl3|) in this case represents a rapidly varying function of x. Therefore the integral in (jl3p 
can be estimated by means of the method of stationary phase. This method can be applied 
for all angles except for a narrow cone of angles, 6k {6k is the angle between the vectors n 
and k) where kL6k < 1- In this region the integral can be considered as constant which 
results in the following form for the energy distribution (before angle averaging!): 

^i(k) ^ ^((rO\ 0k<6o^ {kL)'\ (14) 



where n is the density of discontinuities (= N/ S) and 



_ 

Tl = T{x)dx, I = X2 — Xi, (1) = L. 

J Xl 



For angles 6^ lying far from = i^L)^^ the integral in (fT^ can be estimated by means of 
the method of stationary phase. However, the leading order, proportional to (kr)^^, gives 
zero input because r(xi,2) = 0. Therefore one needs to keep the next order approximation 
that gives: 

- i^ (k.nm.rr » 

where T' is the derivative of F taken at the endpoints Xj. This formula demonstrates singular 
behavior for e(k) at angles 6^ close to and 7r/2 (as well as, to vr and —7t/2). At small angles 
6k < (kL)'^ this expression has to be matched with (fT^ . For the angle range close to 7r/2 
the integral in (fT^ should be cut-off due to the bending of the line of discontinuity. This 
factor switches on at angles — 7r/2| ~ {ka)~^ where a is a characteristic bending length of 
discontinuity (roughly of the order of L). Thus, the energy density distribution e(k) has a 
very narrow angle maximum at 6k near zero with decay at large wave numbers as ~ this 
results in the energy spectrum E{k) ~ k^^, which, up to the logarithmic factor, corresponds 
to the Kraichnan spectrum . For all other angles e{k) decays proportionally to k~^ at large 
k. 

We would like to stress once more that the formulas (|14j) and (jisp are the results of 
non-complete average, i.e, the average with respect to coordinates xia,X2a,yoa- In order to 
get the final answer for the energy spectrum it is necessary to average with respect to angles. 
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E,{k) = 2k e,{k)dd, ^ ^-^((n)'). (16) 



Let us assume first that the angle distribution is isotropic. Then, integrating over angles 
it is easily seen that from the first region (fT^ we have the following contribution: 

'do 

Here the factor 2 appears because of two equal contributions from two regions near 6^ = 
and 6k = vr. Averaging ()15|) over angles gives divergence at — and 6 7i/2. The main 
contribution to the energy spectrum comes from the cut-off at small angles ~ {kL)~^\ 

Thus, both regions of angles give contributions of the same order of magnitude. The com- 
plete answer for the energy spectrum for isotropic turbulence (i.e., isotropic distribution of 
discontinuities) is given by the sum of (fTBj) and (fT7|): 



E{k) 



2Ti^k*L 



((n)Vf-((n^) 



(18) 



which coincides with the spectrum obtained by Saffman 

In order to find the spectrum in the anisotropic situation one needs to average expressions 
()14|1 . (fTH|l with the corresponding distribution function. In numerical experiments anisotropy 
can be conditioned by box boundaries as well as by anisotropy of the pumping of turbulence. 
In the case when such ordering is strong enough the spectrum may get some peculiarities 
originating from non-averaged spectra given by (fT^ . ()15|) . If the width of the angle distri- 
bution function will be narrower than 9q = (kL)^^, then in the angle cone 6^ < A6 the 
energy spectrum E{k, 9) will have the fall-off ~ A;~^, i.e., the same power dependence as for 
the Kraichnan spectrum. Note, however, that this asymptotics is only intermediate because 
^0 = {kL)~^ decreases with increasing k. Therefore starting from k = k* = {L6q)~^, the 
average over angles becomes sensitive relative to the singularities of (fT^ that results in the 
spectrum decreasing proportional to the Saffman fall-off. If the influence of anisotropy is not 
so essential then we should expect the spectrum close to the Saffman one, of course, in the 
case when the main contribution to the spectrum at large k is connected with discontinu- 
ities. The most interesting observation following from the analytical results of this Section 
is that in the very anisotropic case with strong ordering of discontinuities the sharp angular 
maximum of the spectrum has the power fall-off coinciding with that for the Kraichnan 
spectrum. While in the isotropic case our answer coincides with the Saffman answer. In the 
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next Section we present the results of numerical simulation of decaying 2D turbulence at 
high-Reynolds numbers. In particular, the appearance of the power law tails in the energy 
spectrum at large k can be explained rather by discontinuities, than by a cascading process 
with constant enstrophy. 



IV. NUMERICAL INVESTIGATIONS 



To support the arguments of the previous sections and reveal the direct connection be- 
tween the formation of the sharp vorticity gradients and the tail of the energy spectrum 
we have performed a numerical study of the evolution of decaying 2D turbulence. The 2D 
Euler equations ([Q) in the vorticity-streamfunction formulation are integrated numerically 
on a double periodic domain by employing a high resolution fully de-aliased spectral scheme: 

|^ + {cu,7/.}=/i2„V'"u;, (19) 

where ip is the streamfunction related to the vorticity by the Poisson equation: u = —V'^ip, 

the velocity is given as v = (t>^, Vy) = Vip x z and the bracket 

. _ ^ du dijj dip duo 
' dx dy dx dy 

In solving we have included a hyperviscosity term on the right hand side of the equation 

to keep the integration scheme stable (typically we have used n = 3 and /ig = 10~^°). This 

term was observed to decrease the energy by less than 0.002% and the enstrophy by less 

than 20%. We verified that the global features of our results were not dependent on the type 

of viscosity (alternatively we used kinematic viscosity). In the present context we apply the 

hyperviscosity to allow an as wide a dynamical range as possible with the given resolution. 

The domain size is taken to be unity and the resolution is 2048 x 2048 modes. For the time 

integration we employ a third order stiffly-stable scheme. We have chosen the time scale to 

correspond to uJq'^, where corresponds to the maximum vorticity. 

As initial condition we have placed a number of positive and negative vortices randomly 

on the domain, ensuring that the total circulation is zero. Vortices of various shapes/profiles 

from vortex patches (Rankine vortices) to Gaussian vortices form the initial condition. In 

the simulation run described here, we have thus used 10 vortices of each sign with Gaussian 

profiles: 

cu{r,e) = ±uJoexp{-ryRl), (20) 
12 



where tuo = 1 for all vortices, while their radii Rq are varying in the range 0.02 < Rq < 0.075. 
In Fig. la we show the initial vorticity field and Fig. lb the vorticity field at T = 95, which 
corresponds to around 8 vortex internal turnover times {T^ = Atc/ujq). The vorticity field 
has the typical structure for 2D turbulence; it consist of large scale structures (vortices) 
with concentrated vorticity and strongly filamented structures between the vortices. At this 
time there is still strong dynamics in the fiow evolution, with shearing and straining due to 
mutual interactions of nearby structures. Corresponding to the vorticity field we show the 
instantaneous one- dimensional energy spectrum E{k) (total energy: E = E{k)dk) in Fig. 
El The spectrum E{k) for T = in Fig. El shows the spectrum of superimposed Gaussian 
vortices. The spectrum is expanding to the high k-values and at T = 95 a k""" spectrum 
is developed at high wavenumbers for the present case a ~ 3, as clearly demonstrated in 
Fig. El which show the compensated spectrum k^E{K) being constant over almost 2 decades 
in k. 

To investigate the details of the dynamics and how the k~°' spectrum is generated, we 
plot in Fig. ini the di-vorticity field B defined in Sec. II, showing the length |B|, which 
is equal to |Vc(j|. It is clear from the figure that very sharp vorticity gradients appears. 
These are localized in stripes that are mostly along straight lines. The stripes are mainly 
formed outside the dominating vortex structures, and their formation can be explained by 
the analysis discussed in Sec. II. Furthermore it is evident that the concentration of the 
stripes are relatively low. Following the time evolution of maximum value of |B|, Bmax, we 
observe a very rapid growth to values more than 100 times the initial value. Bmax oscillates 
in time with a typical period related to the vortex turnover time, T^. The highest maximum 
reached during this simulation approach 1000, which with a maximum value of the vorticity 
ujo = 1 corresponds to the width of the filaments 6 < 0.001, it is evident that the growth of 
|B| is arrested by the hyperviscosity, and indeed Bmax scales with ij,q^^^. We compare the 
structure of the di-vorticity field with the high pass filtered vorticity field shown in Fig. 4. 
The very similar structure of the high pass filtered vorticity field and the vorticity gradient 
field strongly suggests that the vorticity gradients are responsible for the large-fc part of the 
spectrum, i.e., the part. 

To further discuss the dynamics we show the Weiss field jl3] in Fig. 5, defined hj W = 
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— uj"^) where 

S = [(^ - + + ^^2^1/2 

dx dy dx dy 

is the rate of deformation. Comparing Figs. 3 and 5 we observe that the vorticity gradient 
stripes are ahgned with the contours of W in the strain dominated regions, > 0, mainly 
at the edge of the vortex structures and in between the structure. A careful inspection, 
however, reveals that in the stripes W < 0, i.e., vorticity dominates. This is in line with the 
original arguments of Weiss (see also jfol, 13) ^^at vorticity gradients will tend to concen- 
trate in the strain dominated regions. In particular in the work of Chen et al. loj it was 
demonstrated that the dynamics leading to the enstrophy cascade in driven 2D turbulence 
is most significant in strain dominated regions. 

The spectra we have observed is characterized by having the exponent close to a = 3. 
Thus, with reference to Sec. Ill this will correspond to the spectrum in the anisotropic regime 
where the stripes of vorticity gradients are near straight lines. Indeed in Fig. 3 we see that 
we have stripes that are close to straight lines and the observed spectrum is this in keeping 
with the expectations. To illustrate the anisotropic nature of the spectrum directly we plot 
in Fig. 6 the two dimensional spectrum, e{kx,ky). We observe a clear anisotropy, which 
become particular apparent in the compensated spectrum in Fig. 6b, which is obtained by 
subtracting the angle average of e{kx, ky) (i.e., {2tt)^^ J^'^ ^{kx-, ky)dO). Here we observe a clear 
angular structure. We should emphasize that the spectra obtained here are instantaneous 
spectra obtained a one time and for one realization. Ensemble averaged spectra are expected 
to show a much higher degree of - if not complete - angle isotropy. 



V. CONCLUDING REMARKS 



We have performed a detailed investigation of the relation between turbulent spectra and 
possible singularities in 2D turbulent flows. First, we have demonstrated that the fc-behavior 
of the spectra generated by sharp vorticity gradients, based on the compressible advection of 
di-vorticity, depends significantly on the anisotropy of the spectra. If the angular spectrum 
distribution has one or more very sharp peaks then the one-dimensional spectrum has a tail 
falling-off like k~^ at large k, which is resembling the Kraichnan spectrum, derived from 
spectral cascade arguments. In the opposite case of an isotropic smooth angular dependence 
the spectrum has an asymptotic behavior k~^ as for the Saffman spectra. These arguments 
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allow us to suggest interpretation of many numerical experiments where the spectral expo- 
nent varies in the whole interval between —3 and —4. For instance, in the simulations by 
Okhitani the e-regions, because of their geometry, would give the main contributions 
to the isotropic component of the spectrum which explain the Saffman exponents for the 
observed spectrum in Q|. For /i-regions the situation is different: these regions contain 
stretched vorticity level lines and their contribution to the spectrum should be expected to 
be very anisotropic. This is why these regions produce the /c-behavior as for the Kraichnan 
spectrum. A similar situation takes place in our simulations in comparison with the nu- 
merics performed before [2^|^|. In the both latter simulations the spectra were isotropic 
resulting in spectral exponents like for the Saffman spectrum. In the simulations presented 
in the present paper we have very strong vorticity gradients concentrated on the very narrow 
stripes and therefore the exponent is close to that for the Kraichnan spectrum. Employing a 
filtering of the vorticity field indicates, more or less one-to-one, that the tail of the spectrum 
originates from the sharp vorticity gradients. A strong amplification of di-vorticity of more 
than hundred times is one of the main results of our simulations. A detailed investigation 
of the growth of the di-vorticity maximum is beyond the main scope of this paper and will 
be considered in future works. In conclusion, we stress that this amplification has a natural 
explanation due to compressibility of the mapping ((Tj) providing the transfer from the Eu- 
lerian description to the system of movable curvilinear di-vorticity lines as described in Sec. 
II. 
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Figures 
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FIG. 1: a) Initial vorticity field, b) Vorticity field at time 95 corresponding to 10 vortex turnover 
times. Red color designates positive vorticity and blue color negative vorticity; maximum value is 
1 and minimum value is -1. = 1 
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FIG. 2: Compensated energy spectrum at different times k^E{k) corresponding to the vorticity 
field shown in Fig. ^ 
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FIG. 3: The di-vorticity field B at T = 95. Here the length of di-vorticity vector with the maximum 
(red) value being 673. 
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FIG. 4: High pass filtered vorticity field from Fig. lb, k > 10. 
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FIG. 5: Weiss field for tlie vorticity field show in Fig. lb. Red designates positive values, i.e., 
strain dominated regimes. Blue designates negative values, i.e., vorticity dominated regimes. 
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